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Abstract 

We construct numerically gravitational duals of theories deformed by localized 
Dirac delta sources for scalar operators both at zero and at finite temperature. 
We find that requiring that the backreacted geometry preserves the original scale 
invariance of the source uniquely determines the potential for the scalar field to 
be the one found in a certain Kaluza-Klein compactification of IID supergravity. 
This result is obtained using an efficient perturbative expansion of the backreacted 
background at zero temperature and is conhrmed by a direct numerical computation. 
Numerical solutions at finite temperatures are obtained and a detailed discussion 
of the numerical approach to the treatment of the Dirac delta sources is presented. 
The physics of defect conhgurations is illustrated with a calculation of entanglement 
entropy. 
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1 INTRODUCTION 


1 Introduction 

Many phenomena in condensed matter physics involve strongly interacting systems 
and it is suspected that in many cases the physics is governed by a quantum critical 
point, with effective scale invariance. Thus the description involves strongly coupled 
conformal field theories (CFT). A widely used tool to analyze certain strongly coupled 
CFT systems is the anti-de Sitter conformal field theory (AdS/CFT) correspondence 
mEiE], a setup in which a nongravitational system is mapped to a theory with 
gravity in a higher dimensional spacetime. Over the past few years, it has been 
applied to model many interesting phenomena, including superconductivity and 
superfluidity, Fermi surfaces and non-Fermi liquids - see HIS] for a review. 

The initial applications of AdS/CFT to model condensed matter like systems 
typically involved exact translation invariance which then by construction made 
some relevant physics inaccessible or obscured. Hence one of the main directions in 
current research is to incorporate such key elements of condensed matter systems as 
atomic lattices and localized defects into a dual gravitational AdS/CFT description. 

There have been various ways of how to introduce lattices into the correspondence 
including Q-lattices HEIEI, helical lattices pun], single momentum modes in scalar 
field (neutral lattices) [HI |T2] or in chemical potential (ionic lattices) [T^ flTl IT5] . 

Also very recent studies of localized charge defects M ini brought interesting 
physical insights an example being holographic Friedel charge oscillations [T6|. Most 
of these models incorporated the breaking of translation invariance through the 
introduction of spatially varying smooth sources deforming the CFT action typically 
including one to a few Fourier modes. 

A quintessential model of a solid state lattice is the venerable Kronig-Penney 
model [T8|, in which a periodic lattice of Dirac delta functions is imposed in the 
potential in the Schrodinger equation. This allows for quasi-analytical calculations of 
the key physical quantities of interest. Our main motivation is to develop techniques 
for dealing with a similar setup on the dual holographic gravity side. Although 
seemingly this looks like an innocuous generalization from a single Fourier mode 
source to one where all the Fourier modes are turned on in the soured in a uniform 
^To avoid confusion, throughout this paper by ‘a source’ we always mean the deformation of the field 
theory CFT Lagrangian. Within the AdS/CFT dictionary this is encoded in the boundary conditions for 
the bulk gravitational/scalar/gauge fields. The bulk Einstein-matter equations do not have in contrast 
any external sources. 
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manner, the nonlinear nature of the bulk gravitational description makes this problem 
very challenging. This is especially so as we need to employ numerical relativity 
methods and imposing distributional Dirac delta like boundary conditions is very 
nontrivial and we lack guidance from the conventional numerical relativity literature. 
By itself this problem is thus also quite interesting purely from the numerical 
relativity point of view. 

This paper is a first step in this direction where we construct a fully backreacted 
gravity and scalar configuration with a single (ID) Dirac delta source both at zero 
and at nonzero temperature. We will consider a periodic version in a forthcoming 
work. 

We should note, however, that investigations of localized configurations in hologra¬ 
phy have, despite recent applications, a much longer history starting with a seminal 
Janus configuration of Type IIB supergravity m- Subsequent research include 
configurations with different supersymetry breaking patterns [201 |2ll [221 Ea [Ml |25] 
and configurations at finite temperature [26] . A particulary interesting case has been 
found in the context of eleven dimensional supergravity m, where localized defect 
solution was demonstrated for a scalar field of mass m? = —2. This solution is a 
one parameter, regular deformation of the AdS^ x 5^ vacuum and preserves half 
of the original supersymmetry as well as residual conformal symmetry. However, 
it contains non trivial profiles of various p—forms and the special ansatz makes it 
difficult for finite temperature and finite charge generalizations. 

This motivates the search for cases more suited for applications. In particular 
we will employ a relatively minimal framework developing on ideas initiated in m 
where a single Fourier mode source for the scalar field played a role of an atomic 
lattice. We consider the opposite limit of exciting all modes with equal amplitudes. 
This means imposing a Dirac delta function source for the scalar field. Surprisingly, 
we find that a scale invariant Dirac delta source (which occurs in our setup for 
a linear ID Dirac delta source) leads to a consistent scale invariant backreacted 
background only if the self-interaction potential for the scalar field is exactly equal 
to the one appearing in dimensional reductions of supergravity. In this case we 
can control the numerics directly for the case of a Dirac delta source without any 
regularization and extend the computation to nonzero temperature. 

The paper is organized as follows. In section]^ we formulate the problem and 
present a linearized analysis. Section contains a construction of the backreacted 
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solution at T = 0 together with the unique determination of the scalar potential 
from the requirement that the backreacted geometry preserves the scale invariance 
of the linear Dirac delta deformation. The finite temperature solution is constructed 
in section using the DeTurck method adapted to the Dirac delta asymptotics. We 
also perform numerical cross checks with regularized delta-like sources. The physics 
of the configuration is illustrated with a calculation of the entanglement entropy in 
section both for the T = 0 and nonzero temperature cases. We close the paper 
with conclusions and a summary. 


2 Defect source for the scalar field 


Let us consider a general action for a real, self-interacting scalar field coupled to 
gravity 


S 


1 


IGttGn 


J 


R - 2 Va0V> - ViR 


( 1 ) 


where we choose 

oo 

v{4>) = -6 - 02 _ ^ ^ (2) 

k=l 

with the coefficients Ck being for the moment arbitrary. With such a definition the 
leading order terms give the right cosmological constant and determine the mass of 
the scalar to be m? = — 2. The vacuum solution of this theory is empty AdS^ space 
with a vanishing scalar. Throughout this paper we consider the Poincare patch 


ds^ 



dR -|- dx'^ + dy^ 


( 3 ) 


with z being the bulk coordinate. The asymptotics of the scalar field near the 
conformal boundary {z = 0) are 


0(X, Z) ~ 0l(x)z + (l>2{x)z^ . 


( 4 ) 


According to the standard holographic dictionary m? = A(A — 3), and hence we 
have two allowed solutions A = 2 and A = 1. Both of the choices are possible 
and we choose to set 0i(x) as a source of an operator 0{x) of conformal dimension 
A = 2. Then the subleading term is related to the corresponding expectation value 
02(a:) = (0(x)). 
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2.1 Linearized analysis 

As a first step of the analysis it is instructive to impose a ID Dirac delta source 


(j)i{t,x,y) = Aod{x) , 


( 5 ) 


and find a linearized scalar profile around empty AdS^ with neglected backreaction 
on the geometry. The solution which has the correct boundary conditions for a delta 
function located on a line x = 0 is 




AqZ^ 

7r(x2 + ’ 


( 6 ) 


which is essentially just a bulk-boundary propagator. On the dual theory side this 
corresponds to the shift of the original lagrangian 


S = cScFTs + 


d^x 4>i{x) 0{x) , 


( 7 ) 


which induces the vacuum expectation value {0{x)) ~ 1/x^. This approximation is 
valid for Aq << 1. 

This deformation plainly breaks the translational invariance of the theory. The 
original SO{3, 2) conformal symmetry is broken to the SO{2, 2) conformal symmetry 
of d = (1 + 1) dimensions along the defect, where the operator is sourced. This is 
easily seen when one remembers that ID Dirac delta has scaling dimension equal 
to 1, which, together with scaling of the sourced operator, exactly cancels scaling 
dimension of the integration measure. In order to make this symmetry manifest 
we can adopt the AdS^ slicing coordinates [iniETj in which the background metric 
takes the following form 


ds^ = 


1 


2 , dr^-dU + dy'^ 


dot T 


( 8 ) 


cos(a)2 , , , 

In these coordinates linearized fluctuations around ([^ are governed by the equation 

2 


d dlin I r) , / \ d4>lin , “ j n 

+ 2 tan(Q:)—-h 77777^17(pUn = 0 , 


(9) 


da cos(a)2 

with the solution possessing the right boundary condition being of the simple shape 




Aq 

TT 


cos(a)^ . 


( 10 ) 


It is easy to see that the coordinate change 2 : = rcos(a) and x = rsin(a) this 
solution transforms back to ordinary Poincare coordinates Q. 
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3 Backreacted solution at T = 0 


In this section we will construct a fully backreacted solution for the one dimensional 
delta-like defect ([^ at zero temperature T = 0. We will require the solution to 
possess the residual scaling symmetry of the linearized case. In order to solve the 
full set of Einstein-scalar equation we adopt the AdS^ slicing coordinates familiar 
from the consideration of Janus solutions with a minor modification 

1 


ds = 


A{a) 


da^ dr"^ — dU + dy'^ 

p2 


( 11 ) 


Here p is a constant introduced in order for the AdS boundary to be located always 
at a = 7r/2. This is ensured by imposing A[tt/2) = 0 as a boundary condition. In 
this coordinate system AdS space with linearized scalar profile Q is 


A{a) = cos(a) , 


P=I , 0iin(a) = — cos(a)^ 

TT 


( 12 ) 


The conformal boundary consists of two parts oq = ±7r/2 joined together along 
the defect. Transformation to the Poincare coordinates is in that case 2 = rcos(a) 
and X = r sin(Q!) and will get modified in the full backreacted solution. Due to the 
symmetry of the problem all the relevant fields depend only on a. The coupled set 
of Einstein-scalar equations reads 


^ah - ^ (Va(/)Vfe(?i-bfl'afeH((/>)) = 0 , 


which in our case explicitly gives 


V (<^(a)) — (^(«) ^A(a)(/)'(a)^ — QA"{a)^ + 6W(a)^^ = 0 , 

A{a)A"{a) — — 2A{a)‘^ — -E((/>(a)) = 0 , 


(13) 

(14) 

(15) 

(16) 
(17) 


p^A{a) {A{a)4>"{a) — 2(j)'{a)A'{a)) — V'{(f){a)) = 0 . 

From the first two of the above equations we can obtain a first order ordinary 
differential equation for the function A{a) convenient for numerical or perturbative 
analysis 

-I- A{aY ^6 - -|- V((/>(a)) = 0 . (18) 

It is instructive to first perform a perturbative expansion of the solution 


"4(«) = ^ An(a)e^” , </>(«) = ^ fn(a)e 


2ri.+l 


P=J^ Pne 


2n 


(19) 


n=0 


n=0 


n=0 
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where the lowest order is the AdS solution with the scalar profile ( 12 ). It is convenient 
to identify the expansion parameter e with the value of the scalar field at a = 0 . 

We demand the following boundary conditions for equations : we assume 

reflection symmetry around x = 0 which implies 


dMO) = 0 , dafniO) = 0 . 


( 20 ) 


For the scalar field in addition we require /n(0) = 0, which ensures that e becomes a 
physical expansion parameter i.e. it remains equal to the value of the scalar field at 
a = 0 at any order in the perturbative expansion. We determine the constants pn 
by requiring that the AdS boundary is always at a = 7r/2 through Anij^jT) = 0 . 

The above conditions, for a given choice of the scalar potential (§ determine 
a unique solution. A surprising generic feature of the obtained solution is that 
fn{^) 7 ^ 0, which, when translated to standard Fefferman-Graham coordinates leads 
to a nonvanishing nonnormalizable mode away from x = 0 i.e. a modification of the 
original Dirac delta source to 

(/>i(x) = e(5(x) + + .. .)t^ . ( 21 ) 

\x\ 

Since we want to have a purely localized Dirac delta source, we impose the additional 
condition 

/; (i) = 0 ’ (22) 

as an equation for the coefficients of the scalar potential Q which turn out to be 
uniquely fixed order by order. The first couple of coefficients are ci = 1/36, C 2 = 
1/3240, C 3 = 1/544320, C 4 = 1/146966400, C 5 = 1/58198694400. Those are exactly 
the same as the first terms in the Taylor series expansion of 


D((/) = - 6 cosh(^A) . (23) 

Consequently, only for this potential does a backreacted conformal defect with a 
Dirac delta function source on a line exist. 

It is interesting to note that the above potential is not accidental and comes from 
a certain Kaluza-Klein compactification of D = 11 supergravity after truncation of 
equations of motion to AA = 2 supersymmetry [28l[29] (see [30] for a review). The 
minimal lagrangian of such a reduction, apart from the scalar and the graviton, 
contains one U{1) gauge field coupled in a non-minimal way with the scalar field 


S 


1 


IOttGtv 


J d^x^/^ 


R - 


1 

2 


Va</>V“(/ + 6 cosh 
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Figure 1: Metric and scalar field for 4‘{0) = 1.2. Numerical solution (points) with 
N = A7 spectral grid. Lines correspond to fourth order perturbative solution. 



The result that for the conformal defect to exist the scalar self interaction 
potential has to be of the SUGRA form is indeed quite surprising at first sight. This 
may be qualitatively understood in analogy to deforming a CFT by marginal or 
exactly marginal operators. In the former case, on the linearized level one still has a 
CFT, but if the operator is not exactly marginal, at higher orders in the deformation 
parameter a mass scale is generated and the deformed theory looses scale invariance. 
The phenomenon that we are seeing here is analogous but for a linear ID Dirac 
Delta source. On the linearized level we have scale invariance just by dimensional 
analysis. However we may expect anomalous scaling on the fully nonlinear level w.r.t. 
the deformation parameter. The fact that this does not happen for a theory with 
a supergravity dual is quite natural as we may expect that for a supersymmetric 
field theory there may be appropriate cancellations which would ensure the ‘exact 
marginality’ of the defect deformation. However it would be very interesting to 
understand this in more detail. 

As a cross check of the perturbative considerations we implemented the system 
([^-([Tt]) numerically for the specific supergravity choice of potential ( [23l ) and checked 
that we can obtain a consistent backreacted geometry with the purely localized Dirac 
delta source (i.e. satisfying = 0) for finite values of (/>(0). It is at this stage 

that the introduction of the constant p was particularly useful as it made the size of 
the numerical grid to be fixed and the same irrespective of the value of </>(0). 

For numerical simulations we used standard spectral collocation method m with 
Chebyshev polynomials to account for the a dependence and solving the resulting 
non-linear algebraic equations by the Newton-Raphson method. Curves on the plots 
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in Fig. [^were calculated with N = 47 spectral points. 


4 The finite temperature solution 

In view of possible applications it is interesting and natural to generalize the 
configurations from the previous section to the finite temperature case. This will no 
longer reduce to a system of ordinary differential equations and will depend on both 
variables a and z making the problem much more involved. We will employ the 
DeTurck method with appropriate modifications for incorporating the Dirac delta 
source. 


4.1 The DeTurck method 


The DeTurck method amounts to adding to the original Einstein equations carefully 
chosen terms that make them elliptic partial differential equations. One solves the 
resulting equations numerically and then makes sure that this solution solves the 
initial problem with the original Einstein equations. It was first used to prove the 
short time existence of solutions of the Ricci flow equation [32]. In the context of 
finding numerically static black hole solutions it was used in references |33t l3l] . 

In the general case the equations of motion take the following form 


Gab = -Safe - - ( Na(l)Nb4> + V((/>)gab j = ^ ’ 


(25) 


VaVV 



(26) 


Gab - V = 0 


(27) 


where 

r = (28) 


and r(g) are Christoffel symbols of the reference metric, which we take to be the 
standard black hole metric. 

The generic ansatz for the metric we take is 


ds = 


+ Si(x, z)(dx + F(x, z)dz)^ + S 2 (x, z)dy‘^ 


(29) 
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where we have set AdS radius to L = 1 and factored out 

G(z) = 1 + z + . (30) 

With proper regularity conditions at z = 1 this geometry will have a smooth horizon 
with the temperature T = This ansatz was used in the study of |llj 

where the boundary source was a single Fourier mode 4>i{x) oc coskx. For the case 
of a Dirac delta source 4>i{x) oc 5{x), we need to judiciously modify the coordinate 
system in order to take into account the high variability of the metric coefficients 
close to Dirac delta source at the boundary x = z = 0. Note that in the DeTurck 
method, since the ansatz for the metric is always the most general, the change of 
the coordinate system essentially amounts to an appropriate modification of the 
reference metric. 

Motivated by the treatment of the T = 0 case we define 

tan(a) = — . (31) 

keeping the second relevant coordinate z unmodified. In contrast to the T = 0 
case the geometry will depend on both variables (z, a) and at z = 1 we assume the 
appearance of a regular, static event horizon. 

4.2 The boundary ODE system 

In the new coordinate system the two sides of the boundary on both sides of the 
defect are represented by two points z = 0 and a = L'kI2. The unknown functions 
on an open interval {0} x (—7r/2,7r/2) represent the (backreacted) infinitesimal 
neighbourhood of the delta source and thus have to be determined from the equations 
of motion. The resulting solution will then provide the right boundary conditions 
at z = 0 for subsequently solving the DeTurck equations in the bulk. This is the 
major necessary modification of the standard setup, where typically the boundary 
conditions at z = 0 are completely trivial and explicitly known from the outset. 

To solve the above problem we expand the equations of motion near the z = 0 
point and take the leading order terms, which provide a closed self-consisten10 set 
of coupled, second order ordinary differential equations for the boundary values of 
the fields which we then numerically solve. For a = 0 we set symmetric boundary 
^In particular no z derivatives appear. 


10 







4.3 The hackreacted geometry for T ^ THE FINITE TEMPERATURE SOLUTION 


condition i.e. we set all functions, except for F(a,0), to be symmetric. The off- 
diagonal function has to be clearly anti-symmetric. For a = 'kI2 we choose empty 
AdS boundary condition with a vanishing scalar. The obtained solution at 2 : = 0 
will be used as a Dirichlet boundary condition at the conformal boundary for the set 
of equations in the bulk. 

4.3 The backreacted geometry for T^O 

As was discussed in the previous subsection we assume symmetry in a variable and 
use solution of the boundary ODE system as a Dirichlet condition for z = 0. The 
line a = tt/2 corresponds to the line x = 00 . We set there AdS black hole condition 
with a vanishing scalar field. At the horizon z = 1 we impose the regularity of the 
solution following m- As in the previous case of zero temperature we use the same 
spectral methods m combined with a scipy nonlinear solver in Python. Sample 
solutions for {a, z) G [0,7r/2] x [0,1] are shown in the plot on Fig. We see that the 
special choice of coordinate system rendered the Dirac delta source well behaved with 
the solution being regular and showing no anomalies which are present in the case 
without supersymmetric potential. The geometry smoothly interpolates between 
the horizon and the delta defect located at the boundary. Surprisingly there is only 
little dependence on the holographic direction as one goes from the UV to IR where 
the horizon is located. 




Figure 2: Metric and scalar field for c/){h) = 1-0. Numerical solution (points) with 
Na = Nz = 35 spectral grid. 


While solving (27) numerically we have to make sure that the solution is also a 
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solution to the original Einstein equations. In the pure gravity case one can show 
that it is indeed true and no Ricci solitons are possible [33]. In the case with matter 
this question remains open and we checked numerically, that our solutions indeed 
have ^ = 0. 


4.4 Regularized Dirac delta at T^O 


To cross-check our numerics we also constructed backreacted geometries correspond¬ 
ing to a set of regularized localized defects, which converge to an exact Dirac delta 
for a vanishing regulator. In this case we take as a quasi-localized source for the 
scalar 

Mx) = , (32) 

which is motivated by the linear solution of our problem. Here zq determines the 
width of the conhguration and in the limit zq 0 we recover the strict Dirac delta 
source of the previous sections. 

In this case we use a modified version of our coordinates namely we define 

tans =- , (33) 

z + zo ’ 

so that our source term takes the form 


0i(s) = 


1 


(34) 


2 ;o(tan^ s -|- 1) 

As previously s is in the range s G [0,7r/2]. 

The treatment of boundary conditions is similar to the previous case but for the 
signihcant simplification that for the z = 0 conformal boundary we can now directly 


impose the source (34) for the scalar and empty AdS metric for the remaining 
functions. 

In order to compare the regularized solutions with the backreacted geometries 
from the previous sections we need to match the relevant parameters which is not 
completely trivial. The solutions of the exact defect are parameterized by the value 
of the scalar field at the point s = 0 and z = 0. One recovers the amplitude by an 
integration of the source term over the conformal boundary 

Z -h Zo 


Aq = lim 


z^O 


f dx — (h{x,z) = lim [ da d)(a, z) , 

J-oo TTZ ^ z^oJ_^j 2 7rzcos(a)^ ' 


(35) 


where the last equality comes from the coordinate transformation (33) and the I/vr 
term comes from normalization. This function turns out to have a hxed point for 
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Figure 3: Metric and scalar field for 4‘{0) = 1.0. Numerical solution (points) with 
Ns = Nz = 35 spectral grid. The checkered surface is the solution for the exact Dirac 
delta, while the remaining three surfaces correspond to regularized defects with the 
regulator being zq = 0.3,0.15,0.025. 


4>{0, 0) = 1 i.e. for that value of the field = 1- In Fig. we see that for smaller 
and smaller values of the regulator, the results converge to the solution corresponding 
to the exact Dirac delta defect. 


5 Entanglement entropy 

Entanglement entropy is an interesting probe of various physical systems, for which 
there exists a clear holographic prescription due to Ryu and Takayanagi [35l 136] 
(see |37! for a review). Apart from being of direct physical interest, we will employ 
entanglement entropy also as an important nontrivial cross-check of our numerical 
relativity constructions of the backreacted geometries with Dirac delta sources. 
Namely, we will link the small size limit of the T 7 ^ 0 entanglement entropy with 
the T = 0 one. Due to the fact that the relevant T 7 ^ 0 and T = 0 background 
geometries were obtained using quite different methods in qualitatively different 
coordinate systems, the agreement of these observables will be a nontrivial check of 
our determination of these backreacted geometries. 

According to [35l |36] , in a holographic setup the prescription to compute en¬ 
tanglement entropy of some region of boundary CFT is the following: one takes 
a closed curve inside of which the region of interest lies and then computes the 
extremal surface (with respect to the bulk metric) whose boundary is the given 
curve. Then the entanglement entropy is given by the area of that surface. One 
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must proceed with some caution and renormalise the observable, as the conformal 
boundary lies ’at infinity’ of the bulk space-time and therefore any surface reaching 
conformal boundary has an inhnite area. In the case of various kinds of defects 
similar calculations of entropy in CFT have already been done |26[ [38l 139] . In 
particular in reference [39] entanglement entropy has been obtained for a number of 
known defect configurations. 

Since our problem admits translational symmetry along the y axis we will only 
be interested in entanglement entropy contained in a strip of boundary theory which 
has the defect in its center. The first issue one encounters in this calculation is how 
are curves that reach conformal boundary at some given values of the x variable 
mapped under the coordinate transformation connecting Poincare coordinates ([^ 
and our slicing ones Q for the T = 0 case, or the transformation given by (31) 
for the black hole case. The two cases need a separate treatment as the employed 
coordinate systems are different. 


5.1 The T = 0 case 

In the T = 0 case we use the generalized angular coordinate system (r, a) with a 
generic rule of transformation to ordinary Fefferman-Graham (FG) patch coordi¬ 
nates z = rf{a), X = rg{a). The functions f{a) and g{a) are no longer simple 
trigonometric functions in the backreacted case which poses some complications in 
the following. 

The problem of finding the relevant extremal surface reduces to determining the 
function r{a) which solves the equations of motion following from extremalization of 
the Nambu-Goto action. The entanglement entropy is then determined by evaluation 
of the action on the solution. With our assumptions NG action takes the following 
form 

Tng = -- / ^2A( s,2V ^(«)^ > (36) 

pr{aYA{aY * 

with the resulting equation of motion 

A(a)r(a)‘^ [p^r''{a) + r{a)j — 2p^A'{a)r'{a) -f r(a)^^ = 0 . (37) 

Since the integral in the Nambu-Goto (NG) action is clearly UV divergent, one 
has to perform a cut-off procedure. In order to combine the results both with and 
without the defect, the regularization has to be done covariantly. This is achieved by 
always expressing the asymptotic part of the metric in FG coordinates and setting 
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the regulator z = e. This then translates back to the cutoff Uc in new coordinates 
by the solution of e = r{ac)f{oic). We thus have to determine, to some degree, the 
coordinate transformation function f{a). 

This can be done by writing, in FG coordinates, the most general form of a 
metric which is consistent with the symmetries of the problem 




dx^ 


ds = + C{zlx)^r- + D{zlx) 


dy^ - dP 


(38) 


Employing now the transformation laws ^ = rf(a), x = rg{a) and comparing with 
our metric ansatz ( |11[ ), one obtains 

f («) _ yi-Aja)'^ g'{a) _ A{a) 

/(a) pA{a) ’ g{a) p^l - A{aY ' 

The two integration constants are determined by the requirement that for a ^ 'k 
we get the empty AdS metric which boils down to two conditions 


C(V2) = - lim 


a^7r/2 g'{a)g{a) 


D{'k/2) = lim 

CX.—^'Kj2 


C{a)g{af + f{a) 


= 1 . 


(40) 


(41) 


It is important to note, that coordinates (38) do not cover the whole spacetime. 


since ^(0) > 1. However, for our purposes, it is enough to consider only the near 
boundary region, where formulas work fine. There, the transformation can be found 
perturbatively, in parallel to the expansion of A{a), and it reduces to the ordinary 
polar coordinates in the absence of the defect, in which case it is regular in the whole 
range of a. 

To get the entanglement entropy we now must evaluate the NG action on the 
solution. Due to translational symmetry along the defect we quote all formulas per 
unit length in the y direction. As it was already mentioned there is a UV divergence 
coming from a simple pole of A(a) at a = tt/2. This part of the integral can be 
separated and estimated analitically to be 


5'div — 


P 


(42) 


(7r/2 - ac)L 

where Uc is a cut-off dehned above by e = r{ac)f{ac). To find the limiting behaviour 
we can now expand this equation near the boundary to obtain e = L/'(7r/2)(7r/2 —etc) 
to get the divergent part expressed by a Poincare path regulator 

P/'(7r/2) 


5'div = 


(43) 
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Figure 4: Regular part of the entanglement entropy for the defect conhguration 
(dots) and for case of empty AdS^ (line). Here the defect amplitude corresponds 
to 0(0) = 2. 


It turns out that for all our solutions p/'(7r/2) = 1, which gives the right UV 
divergence. The regular part follows conformal invariance so that the total result 
reads 

=- Y , (44) 

e L 

for some positive B which depends on the amplitude of the source. 

A surprising feature of the obtained result is that the regularized entanglement 
entropy is lower than the corresponding one in empty AdS spacetime (see figure]^. 
Indeed for 0(0) = 1, a ht to the difference between numerical data and the empty 
AdS result yields 

-^defect - *S'AdS = -- • (45) 

This is clearly intriguing, although it does not indicate any contradiction, es¬ 
pecially as we are not considering here a different state in the original theory but 
rather a deformation of the theory by a nontrivial source added to the Lagrangian. 

5.2 The finite temperature case 

The entanglement entropy calculation at nonzero temperature follows the same basic 
steps as before but with appropriate modifications due to the different coordinate 
system employed for the background geometry. In this case the corresponding 
minimal surface will be described by (a, y) coordinates with a non-trivial dependence 
z(a). The boundary condition for the surface is z'{yIT) = —L where L is the half 


16 







5.2 The Gnite temperature case 


5 ENTANGLEMENT ENTROPY 


width of the strip. In order for the equations to be regular we adopt the following 
ansatz 

z{a) = Z(a)cos{a) , (46) 


after which the NG Lagrangian takes the following form 

(a, cos{a)Z{a)) 


Lng = 


cos{a)‘^ Z (a)'^ — cos^{a)Z(a)^ 


Va , 


(47) 


where 


>4 = ^1 — cos^(a)Z(a)^^ 5i(a, cos(a)Z(a)) |^Z'(a)(cos(a)F(Q!, cos(Q!)Z(a)048) 

+ sin(a)) + Z(Q;)(cos(a) — sin(a)F(a,cos(a)Z(a))) I + 

+ H 2 {a,cos{a)Z{a))l sin(Q:)Z(a) — cos(a)Z'(a) 



The equation of motion is a non-linear ordinary differential equation defined numeri¬ 
cally in terms of the prohles calculated in section The boundary conditions for the 
extremal surface translate to: Z'{0) = 0 and Zii^jT) = L. The resulting equation 
with these boundary conditions can be solved using spectral discretization with the 
Newton method. 

In order to evaluate the entanglement entropy we regularize the divergent action, 
evaluated on-shell, by subtracting pointwise the proper AdS NG Lagrangian evalu¬ 
ated on a corresponding extremal surface in empty AdS, found numerically in an 
analogous coordinate system on the same numerical grid in a i.e. 


5'def(-^) = [ da{C]^G - TacIs) • (49) 

Jo 

We calculate the entanglement entropy for a Dirac delta defect characterized by 
(/>(0) = 1 for the (half-)width of the strip the range 0.02 < L < 1.1. The results are 
shown in figure 

Due to the conformal symmetry of the defect the generic form of the entanglement 
entropy has to be 


<S'def(-^") 


B{LT) 

Z 


(50) 


where B[LT) is some smooth function. Of course B{LT) will also depend on the 
(dimensionless) amplitude of the Dirac delta function. This formula allows for making 
a connection with the entanglement entropy computation for T = 0 presented earlier. 
Indeed, for fixed T when L —)• 0 we recover the vacuum defect result. For this 
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Figure 5: Entanglement entropy for a Dirac delta defect with 4>{0) = 1.0 at finite 
temperature with = 50 spectral grid. 


particular case the 1/L falloff has a coefficient a = —0.0107 which is consistently 


checked against zero temperature calculation (45). This is also a nontrivial check 
of our numerical constructions of the relevant backgrounds, as due to the different 
coordinate systems employed the extremal surfaces for small L obtained in both 
cases are quite different as can be seen on Fig. 




Figure 6: Plots of the minimal surfaces embeddings for L = 0.05 for T = 0 (left 
panel) and T > 0 (right panel). 


It is quite interesting to compare the finite temperature entanglement entropy 
with and without the defect. The relevant plots are shown in figure Again we 
observe the rather surprising result that the defect source decreases the entanglement 
entropy. Clearly, for small L the effects of the defect are significant while for L > 0.25 
thermal effects start to dominate. 
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Figure 7: Regular part of the entanglement entropy for the defect conhguration 
(blue dots) and for AdS black hole case (red dots). 

6 Conclusions 

In this paper we constructed holographic duals of strongly coupled three-dimensional 
CFT’s deformed by a localized Dirac delta source. The motivation for this work 
was to move towards a holographic construction mimicking a crystalline lattice with 
pointlike localized sources. In the present paper we have concentrated on developing 
the necessary numerical relativity methods in order to consistently handle Dirac 
delta like sources and considered explicitly a single defect along a line both at zero 
and at nonzero temperature. 

The ID Dirac delta source for a scalar operator of dimension A = 2 is scale 
invariant. We found that requiring that the backreacted geometry respects this 
scale invariance imposes very stringent constraints on the scalar potential which 
is consequently uniquely determined. An intriguing outcome is that the resulting 
potential is exactly the scalar potential appearing in certain Kaluza-Klein reductions 
of IID supergravity. All further considerations in the present paper were performed 
with this concrete choice of the scalar self-interaction. 

In order to hnd the dual backgrounds, we had to use systems of coordinates 
which were adapted to the presence of the Dirac delta source and which took into 
account the high variability of the scalar held and of the metric coefficients in the 
inhnitesimal neighbourhood of the point of insertion of the Dirac delta function on 
the boundary. 

At zero temperature, we constructed the dual backgrounds in two ways: using a 
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perturbative expansion and performing a direct numerical solution of the equations 
of motion using an AdS slicing analogous to the one used for obtaining the Janus 
solution [ISIET]. 

For nonzero temperature, we adopted the DeTurck method which required, 
however, two modifications. Firstly, the adapted choice of coordinates (similar 
but different from the one that we used for T = 0) was encoded in the choice of 
coordinates used for the reference AdS black hole metric. Secondly, the boundary at 
z = 0 of the numerical grid represented really the infinitesimal neighbourhood of 
the Dirac delta source and the values of the fields there had to be determined from 
the equations of motion. With those two modifications in place we constructed the 
numerical background and performed two cross checks. We compared the resulting 
numerical background with analogous geometries obtained for regularized delta 
sources. We also compared the small size limit of entanglement entropy with the 
entanglement entropy evaluated at zero temperature. 

Incidentally, we found that the entanglement entropy evaluated in the theory 
deformed by the Dirac delta source is lower than the analogous quantity without 
the defect. It would be very interesting to understand this property from a more 
physical perspective. 

There are numerous directions for further study like the construction of a lattice 
of Dirac delta defects in order to study phenomena analogous to m , extensions 
to chemical potentials and pointlike sources. It would be interesting to determine 
whether the ‘hovering black hole’ phenomena observed in m have their counterpart 
in the present context. We intend to address at least some of these issues in 
subsequent work. 
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